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Presented in this work is the most comprehensive INTEGRAL AGN sample: it includes 272 
active galaxies for which optical identifications and redshifts are available. 2-10 keV fluxes and 
column density measurements are also collected for all the sample sources. For 33 new hard X- 
ray discovered AGN with no previous X-ray coverage, XRT and XMM data analysis is reported 
for the first time in this work. Examples of future studies which are being developed using this 
large X-ray selected sample are presented together with some early results. 
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1. The INTEGRAL AGN sample 

From 723 sources listed in the 4th INTEGRAL/IBIS survey [1] we took a census of 234 
AGN considering only the extragalactic objects for which a firm optical identification and redshift 
measurement is available. In order to a have a complete view of INTEGRAL extragalactic sky we 
have also cross-correlated the 4th IBIS catalogue by Bird et al. 2009 with the INTEGRAL all-sky 
survey by Krivonos et al. [4] updated on the web site, and among their 211 objects identified as 
AGN, we found 38 extra sources not included in our list. With the addition of these extra objects 
the final INTEGRAL AGN catalogue presented in this work lists 272 galaxies (last update March 
2011). Note that new sources are continuously coming in [6]. In figure 1 (left) the whole sample 
is reported in the classical 20-100 keV luminosity vs redshift plot. This sample represents the most 
comprehensive INTEGRAL AGN catalogue and this has two strengths. First: all sources, except 
for one BL LAC object, have optical spectra which means secure identifications and measured 
redshifts; and second: all sources (except for one case) have X-ray data and so 2-10 keV fluxes and 
N# estimates available. 

The 272 AGN belong to the following optical classes: 132 objects are Seyfert of type 1 (Seyfert 

1, 1.2, 1.5 and NLS1) , 116 are Seyfert of type 2 (Seyfert 1.9-2, and XBONG) and 22 are blazars 
(BL Lac and QSO). 

In order to assign to each object the most appropriate optical class, a cross check between the 
class reported in the NED archive and that reported in the last edition of the Veron-Cetty, Veron 
extragalactic catalogue [9] has been performed with the exception of those objects for which the op- 
tical classification has been provided by our own identification works [6, and reference therein]. In 
case of different classification the more recent one has been preferred and for subclass assignments 
(Syl.2, Syl.5, Syl.8, Syl.9) the diagnostic criteria of Winkler [10] have been adopted. 

2. X-ray data analysis and results 

Among the 272 AGN there are 33 sources which are new hard X-ray emitters (discovered by 
INTEGRAL/Swift) and therefore they had no X-ray (2-10 keV) coverage prior to this work. For 
these sources we analysed X-ray data available from Swift and XMM-Newton observations. 

Swift-XRT data are available for 31 sources and they have been analysed in this work for the 
first time. XMM data of RX J0137.7+5814 and IGR J13415+3033 are also considered as, for the 
first source, no X-ray published data have been found, and for the second, a re-analysis has been 
performed getting quite different results form those obtained by [7] using XRT data alone. 

XRT data reduction was accomplished using the XRTDAS standard data pipeline package 
(xrtpipeline v. 0.12.6), in order to produce screened event files. Furthermore, to univocally identify 
the X-ray counterpart, an accurate imaging analysis has been performed for those new AGN first 
detected at high energies (>20 keV) by INTEGRAL/IBIS. For each observation we analysed the 
0.3-10 keV image to search for sources detected at a confidence level > 3a both within the 90% 
and 99% IBIS error circles. If more sources were present within these IBIS uncertainties, we 
restricted the image analysis to the 3-10 keV band to spot those sources with the hardest spectra 
(i.e. those with detection above 3 keV), since these are most likely to be the counterparts of IBIS 
objects. Then, we estimated the X-ray positions using the task xrtcentroid, providing in this way 
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Figure 1: Left panel: Hard X-ray luminosity vs redshift for all the INTEGRAL AGN sample. Circles are 
type 1 objects, squares are type 2 and triangles are blazars. Open symbols are objects where no intrinsic 
absorption have been measured. Right panel: Column density distributions in type 1 (up) and type 2 (down). 
Dashed bins represent upper limit measurements including Galactic values. 

for each source the position at arsec level accuracy in order to be sure of the optical identification. 
Events for spectral analysis were extracted within a circular region of radius 20", centered on the 
source position, which encloses about 90% of the PSF at 1.5 keV, while the background was taken 
from various source-free regions close to the X-ray source of interest. In all cases, the spectra 
were extracted from the corresponding event files using the XSELECT software and binned using 
grppha in an appropriate way, so that the % 2 statistic could be applied. We used version v.011 
of the response matrices and create individual ancillary response files arf using xrtmkarf v. 0.5.8. 
The energy band used for the spectral analysis, performed with XSPEC v. 12.5, depends on the 
statistical quality of the data and typically ranges from 0.3 to 8 keV. 

For two sources of our sample we use X-ray data acquired with XMM-Newton and they were 
reprocessed using the Standard Analysis Software (SAS) version 9.0.0 employing the latest avail- 
able calibration files. Only patterns corresponding to single and double events (PATTERN<4) were 
taken into account and the standard selection filter FLAG=0 was applied. For each observation, we 
analysed the XMM-Newton EPIC (MOS plus pn) images to search for X-ray sources which fall in- 
side the IBIS error box and are therefore likely IBIS counterparts. Next we obtained X-ray spectra 
in the 0.5-12 keV band. Source counts were extracted from circular regions of radius 20" centered 
on the source while background spectra were extracted from circular source-free regions of 80" 
radius. The ancillary response matrices (ARFs) and the detector response matrices (RMFs) were 
generated using the XMM SAS tasks arfgen and rmfgen while the spectral channels were rebinned 
in order to achieve a minimum of 20 counts in each bin. 

In order to estimate the X-ray spectral parameters of these newly discovered sources, we fit the 
X-ray spectrum using a baseline model composed by a primary continuum modeled with a power 
law absorbed by the Galactic column density calculated through XSPEC [2]. To measure the 
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amount of cold absorption intrinsic to the source, an extra column has been added when required. 
In a few cases, due to the poor statistical quality of the data, the photon index was fixed to T = 
1.8 which represents a canonical AGN value so as to allow a measurement of the intrinsic column 
density. In table 1 the data analysis results are reported together with the optical classification of 
each object; the last two are those for which XMM data have been reduced. 

The set of AGN analysed for the first time here, are well representative of the entire IN- 
TEGRAL sample. A few peculiar objects emerge from table 1, like IGR JO 145+6437 which is 
optically classified as Sy2 but shows no intrinsic absorption, and IGR J03249+4041 which is a pair 
of galaxies both at ~0.047, both classified as Sy 2 and both emitting at hard X-ray energies. 

3. INTEGRAL AGN sample: Exploitation Studies 

3.1 X-ray absorption versus Optical classification 

The unified model for Seyfert galaxies explains their different observed characteristics with 
their orientation with respect to the observer line of sight. An AGN is classified as type 1 if there 
is the presence of broad permitted lines in its optical spectrum while it is classified as type 2 if it 
does not show these broad components. In the type 2 Seyfert the broad line region is hidden by 
absorbing material which should largely correspond to that measured in the X-ray band. Therefore 
heavy absorption is expected in type 2 and less or even none in type 1 . It is well known that there is 
a significant number of AGN for which the expected optical and X-ray characteristics do not match 
and we know cases of AGN of type 1 where strong absorption is measured and, on the other hand, 
there are also AGN showing in their optical spectra only narrow emission lines but displaying only 
mild or even absent absorption in their X-ray spectra. 

The INTEGRAL sample of AGN is an ideal tool to investigate this issue since the hard X-ray 
selection is almost unbiased with respect to the absorption. To this end we have performed the 
column density distribution dividing our sample sources in type 1 (including also Syl.2-Syl.5) and 
type 2 (including also Syl.9). The histograms of N# in type 1 and type 2 Seyfert galaxies in the 
INTEGRAL AGN sample is shown in figure 1 (right). If we assume the dividing line between non- 
absorbed and absorbed AGN lies at N# = 10 22 cm~ 2 f 1 ], it is clearly visible from figure 1 (right) 
that a number (23) of type 1 objects show intrinsic absorption higher than 10 22 cm~ 2 , and a number 
(11) of type 2 have column density lower than 10 22 cm~ 2 or even none in excess to the Galactic 
one (dashed bins in figure). If we take into account only the 'pure' Seyfert galaxies (i.e. excluding 
LINERS, XBONG), the total number of AGN which do not fit into the unified scheme are 33/250 
which correspond to 13% of the sample and this percentage is quite close to that typically found in 
the soft X-ray selected samples of AGN [3]. 

3.2 X-ray absorption versus host galaxy axial ratio 

It is becoming increasingly clear that the study of AGN is strongly connected with the prop- 
erties of the host galaxies and, in particular, the inclination of the host galaxy can play a role in 
relation to the X-ray absorption. The inclination is determined by the axial ratio which is defined 
as the semi-minor (b) over semi-major axis (a) so that b/a=l corresponds to a face-on galaxy and 

1 it corresponds to the column density needed to hide the BLR for clouds having a standard gas-to-dust ratio [8] 
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Results of the spectral analysis of the XRT/XMM data 



Source 


Class 


N H (intr) 

(10 22 cm- 2 ) 


r 


F(2-\0keV) 

(\Q- n erg cm- 2 s" 1 ) 


X 2 /v 


IGR J00158+5605 


Syl.5 




2.30 ±0.34 


0.81 


3.7/6 


IGR JO 1545+6437 


Sy2 




[1.8] 


0.0097 


3.8/3 


RBS 345 


Syl 




1.62+024 

—0.22 


0.26 


17.2/8 


IGR J03249+4041-SWW 


Sy2 


> 0.15 


1-53 r\ o n 

0.39 


0.12 


7.56/9 


IGR J03249+4041-NEW 


Sy2 


3 04+ L54 

99 


[1.8] 


0.09 


4.35/6 


IGR J03334+3718 


Syl.5 


0.14+^ 

—0.06 


1.86+0-!? 

—0.17 


0.61 


49.1/35 


LEDA 15023 


Sy2 


30.1+1] 


2.21+0-^ 

—0.50 


0.13 


8.1/8 


LEDA 075258 


Syl 


0.0057+0™ 2 

— D.IK) I 


i 77 +0.09 
—0.09 


0.55 


109.8/109 


IGRJ06058-2755W 


Syl.5 




i 71 +0.09 

—0 09 


0.66 


39.0.50 


IGR 108557+6420^ 


likely Sy2 


19.6+H 

— 4.0 


[1.8] 


0.28 


6.3/8 


IGR J08558+0814 


Syl 




2.24 ±0.46 


0.0065 


4.6/13 


IGR J09253+6929 


Syl.5 


14.8+ 2 f° 

— 1 l.U 


[1.8] 


0.045 


0.22/2 


IGR J09446-2636 


Syl.5 




1.66 ±0.09 


0.50 


38.8/32 


IGRJ13 187+0322 


QSO 





[1.8] 


0.015 


3.0/6 


IGRJ14175-4641 


Sy2 


7 c 5T+283.0 
—40.0 


[1.8] 


0.095 


3.0/5 


IGRJ14301^1158 


Sy2 


1.4+08 

—0.5 


1.83+050 

—0.58 


2.70 


5.5/11 


IGR J1531 1-3737 


Syl 




1.73±0.16 


1.86 


8.4/11 


IGR J 15549-3740 


Sy2 


5.1 +2 -0 

— 1.3 


[1.8] 


0.29 


8.5/9 


IGRJ16426+6536 


NLSy 1 




[1.8] 


0.0085 


4.3/6 


SWIFT J 1650.5+0403 


Sy2 


— 1.3 


, O7+0.21 

1.0/ n "in 

—0.37 


0.32 


23.6/32 


IGRJ17036+3734 


Syl 




2 3O+0- 32 

JU -0 21 


0.21 


4.2/5 


IGRJ17476-2253 


Syl 


— 


1 Q1+1.00 
i-5'i_091 


0.26 


2.7/3 


1RXS J175252.0-053210 


Syl.2 


- 


i -^ J -1.00 


0.26 


2.3/2 


IGR J1831 1-3337 


Sy2 


1 4+5-0 


[1.8] 


0.026 


3.2/6 


IGR J 19077-3925 


Syl.9 


14+°- 10 


, 7R+0.24 
1 - /8 -0.22 


0.36 


11.4/18 


IGR J191 18-1707 


LINER 


10.4^-1 


[1.8] 


0.14 


2.1/4 


PKS 1916-300 


Syl.5/1.8 




2.04 ±0.07 


0.62 


57.7/66 


1RXS J192450.8-291437 


BL Lac 


0.073 ±0.03 


1.95±0.10 


0.81 


104.1/89 


QSOB 1933-400 


QSO 


n+ 013 


1 84+0- 19 

1 '°^-0.31 


0.17 


11.7/13 


IGRJ19491-1035 


Syl.2 


i-^ -0.73 


i 79+0.41 


0.56 


19.68/14 


SWIFT J2044.0+2832 


Syl 




i 94+0.05 

i> ^-0.05 


0.71 


109.1/115 


1RXS J21 1928.4+33359 


Syl.5 


46+ 08 


, 78+O.I8 
1 - /8 -0.17 


0.59 


80.2/61 


RX J0137.7+5814 


BL Lac 




7 , 0+0.04 
^ -0.04 


1.01 


224/210 


IGRJ13415+3033 


Sy2 


29.70+]^ 


i 9 7 +0.i6 


0.25 


426/397 



' interacting galaxies, here a tentative spectral analysis has been performed 
" Best-fit model includes an absorption edge at ^0.52 keV 

b Best-fit model includes a second power law of F= 1.8 (fixed) to account for the excess observed 
below 1 keV. 
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b/a~0 to an edge-on galaxy. We have collected from the literature all the available axial ratios 
for the entire INTEGRAL AGN sample and we have found a different behavior between type 1 
and type 2 objects, with only type 1 AGN showing a deficit at low values of b/a. This has been 
investigated in terms of absorption suggesting the existence of an extra absorbing component on a 
larger scale than the torus assumed in the unified theory. 

3.3 X-ray absorption versus galactic structures and interactions 

Another possible connection between the X-ray absorption and the large scale structures could 
be represented by the presence of large scale gravitational torques such as in barred and interacting 
galaxies as they can transport gas close to the center of the AGN and this gas concentration can 
play a role in the X-ray obscuration [5]. Using our sample we do not find any correlation between 
absorption and the presence of bars and/or interactions, despite the fact that '--40% of our galaxies 
are strongly barred and 38% interacting. 

These are some of the results that can be obtained with our large sample, highlighting the 
potential of statistical and population studies of hard X-ray selected AGN. Further works are in 
progress. 
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